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Abstract 

In a numerical Monte Carlo simulation of SU(2) Yang-Mills theory with 
dynamical gauginos we find evidence for two degenerate ground states at the 
supersymmetry point corresponding to zero gaugino mass. This is consistent 
with the expected pattern of spontaneous discrete chiral symmetry breaking 
Z 4 — > Z 2 caused by gaugino condensation. 
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1 Introduction 



The basic assumption about the non-perturbative dynamics of supersymmetric Yang- Mills 
(SYM) theory is that there is confinement and spontaneous chiral symmetry breaking, 
similar to QCD JT|. (For a more recent introduction and review see also 0.) In the 
past years there has been great progress in the understanding of the non-perturbative 
properties of supersymmetric gauge theories, in particular following the seminal papers 
of Seiberg and Witten ||. In case of N = 1 SYM theory the non-perturbative results 
are not rigorous but fit into a self-consistent plausible picture of low energy dynamics 
of supersymmetric QCD (SQCD) [|J. The features of the low energy dynamics, like 
symmetries and bound state spectra, are formulated in terms of low energy effective 
actions || ||. Lattice Monte Carlo simulations may contribute by directly testing some 
of these predictions. 

The expected pattern of spontaneous chiral symmetry breaking in SYM theories is 
quite interesting: considering for definiteness the gauge group SU(iV c ), the expected sym- 
metry breaking is Z 2 n c — ► Z 2 . This is because the global chiral symmetry of the gaugino 
(a Majorana fermion in the adjoint representation) is anomalous. The symmetry trans- 
formations are 

y x e -^m x , ^ _ y x e-*™ , (1) 

where the Dirac-Majorana fields are used which satisfy, with the charge-conjugation Dirac 
matrix C, 

V x = C* x , % = m T x C . (2) 

The group of symmetry transformations in ([!]) coincide with the i?-symmetry and hence 
will be called U(1) R . The transformation is equivalent to the transformation of the gaugino 
mass rrig and a shift of the ^-parameter: 

m g -> ruge' 2 ^ 5 , 6^6- 2N cV . (3) 

Since 9 is periodic with period 2n, in the supersymmetric case with m g = the U(1) R 
symmetry is unbroken if 

ip = if k = — } (fc = 0,l,...,2jV c -l) . (4) 
Gaugino condensation means a non-zero vacuum expectation value 

= (\ a x \ xa + K\ x& ) ± o . (5) 

(Here, besides the Dirac-Majorana field, the Weyl-Majorana field components A" and A" 
are also introduced.) The gaugino condensate is transformed under U(1) R according to 

^ (% x e- 2i ^ x ) . (6) 
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If such a condensate is produced by the dynamics then it breaks the Z 2 n c symmetry to 
Z 2 : the expected spontaneous chiral symmetry breaking is Z 2 ^ c Z 2 . This implies the 
existence of N c discrete degenerate ground (vacuum) states with different orientations of 
the gaugino condensate according to (f|), @. 

A non-zero gaugino mass (m g 7^ 0) breaks the supersymmetry softly. As a function 
of the gaugino mass the degeneracy of the N c ground states is resolved. At m g = the 
lowest ground state is changing. This gives rise to a characteristic pattern of first order 
phase transitions. 

In the special case of SU(2) gauge group, which will be considered in this paper, we 
have Z 4 — > Z 2 and in the two vacua the gaugino condensate has opposite signs. At m g = 
the lowest ground states are exchanged and a first order phase transition occurs. In this 
letter we report on a large scale numerical Monte Carlo simulation with the aim to find 
numerical evidence for the existence of this phase transition. 

2 Lattice formulation 

The definition of an Euclidean path integral for Majorana fermions may be obtained 
by starting from the well known Wilson formulation || of a Dirac fermion in the adjoint 
representation. If the Grassmanian fermion fields in the adjoint representation are denoted 
by ip r x and ip x , with r being the adjoint representation index, then the fermionic part of 
the lattice action can be written as 



xu,yv 

Here the fermion matrix Q is defined by 

4 

Qyv,xu = Qyv,xu[U~\ = 0~yx&vu K ^ ] (1 T/i)^ou,a;/i &y+fi,xiX T/tt) ^vu,yn ' (®) 

K is the hopping parameter and the matrix for the gauge-field link in the adjoint repre- 
sentation is defined as 

Vrw = V„[U] = 2Tr(U^T r U^T s ) = V r * s> ^ = V~% . (9) 

The generators T r = |A r satisfy the usual normalization Tr (A r A s ) = \5 TS . In case of 
SU(2) we have T r = ~r r with the isospin Pauli-matrices r r . Starting from the Dirac 
fermion fields one can introduce two Dirac-Majorana fields vl/^ 1 ' 2 ) satisfying (fj): 

4rW = ±=y> + cf) , ¥^ = -^(-V, + cf) (10) 
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and Sf can be rewritten as 
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S/ = jEE*S J Ve. (ii) 

Z j=l xu,yv 

Using this, the fermionic path integral for Dirac fermions becomes 

2 

j [##]e~ 5 ' = / [##]e^ w = det Q = f[ J [d^]e~^ Q ^ U) . (12) 

For Majorana fields the path integral involves only [d\E^], either with j — 1 or j — 2 
hence, omitting the index (j), we have 



= J [d^]e-^ T ° QiS = Pf(CQ) = Pf(M) . (13) 



Here the Pfaffian of the antisymmetric matrix M = CQ is introduced. The Pfaffian 
can be defined for a general complex antisymmetric matrix M a p = —Mp a with an even 
number of dimensions (1 < a, (3 < 2N) by a Grassmann integral as 

Pf (M ) = J [d<f>}e-^ M ^ = ^e aiPl ... aN0N M a ^ . . . M Qjv/3iV . (14) 

Here, of course, [d(j)\ = d<p2N ■ ■ ■ d<pi, and e is the totally antisymmetric unit tensor. It can 
be easily shown that 

[Pf(M)] 2 = detM . (15) 

One way to prove this is to use detM = detCQ = detQ and eqs. (p^2|)-(p^). Besides the 
partition function in (|12"D, expectation values for Majorana fermions can also be similarly 
defined 1, ^ . 

It is easy to show [[Hp that the adjoint fermion matrix Q has doubly degenerate real 
eigenvalues, therefore det Q is positive and Pf(M) is real. Omitting the sign of Pf(M) 
one obtains the effective gauge field action [|T^] : 

Scv = /?£ (l - \ttU p1 ) - ~logdetQ[C/] , (16) 

pi 

with the bare gauge coupling given by /3 = 2N c /g 2 . The factor | in front of log det Q tells 
that we effectively have a flavour number Nf = | of adjoint fermions. The omitted sign 
of the Pfaffian can be taken into account in the expectation values: 

(AsignPf(M)) cy 
{ ; (signPf(M)) cy • 1 } 

This sign problem is very similar to the one in QCD with an odd number of quark flavours. 
The value of the Pfaffian, hence its sign, can be numerically determined by calculating 



an appropriate determinant [0] . It turns out that in updating sequences with dynamical 
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gauginos configurations with positive Pfaffian dominate. This is shown by explicit evalu- 
ation on 4 • 8 lattices. It is plausible that the sign chang function of the valence 
hopping parameter, typically occur at higher values than the value of K in the dynamical 
updating [13[]. Therefore, in the present work, we consider the effective gauge action in 
flT6|) and neglect the sign of the Pfaffian. To take into account the sign is possible but 
numerically demanding, therefore we postpone it for future studies. 

Since the Monte Carlo calculations are done on finite lattices, one has to specify 
boundary conditions. In the three spatial directions we take periodic boundary conditions 
both for the gauge field and the gaugino. This implies that in the Hilbert space of states 
the supersymmetry is not broken by the boundary conditions. In the time direction we 
take periodic boundary conditions for bosons and antiperiodic ones for fermions, which is 
obtained if one writes traces in terms of Grassmann integrals. (The minus sign for fermions 
is the usual one associated with closed fermion loops.) Of course, boundary conditions 
do not influence the physical results in large volumes. For instance, we explicitly checked 
that the distribution of the gaugino condensate is not effected if in the time direction 
periodicity is assumed for the fermions, too (see below). Another interesting possibility 
would be to consider twisted boundary conditions jnj which are useful in theoretical 



considerations about supersymmetry breaking [15 



3 Monte Carlo simulation 

The expected first order phase transition at zero gaugino mass should show up as a jump 
in the expectation value of the gaugino condensate ([|). The renormalized gaugino mass 
is obtained from the hopping parameter K as 



Z m {a\i) 



1 1 
K~K~ 



= Z m (ajj)m g . (18) 



Here a denotes the lattice spacing, \i is the renormalization scale and K = Kq{(5) gives the 
/^-dependent position of the phase transition, which is expected to approach K = 1/8 in 
the continuum limit (3 — > oo. The bare gaugino mass m § is defined, as usual, by omitting 
the multiplicative renormalization factor Z m . The renormalized gaugino condensate is 
also obtained by additive and multiplicative renormalizations: 

(^x^ x )r M = Z(afi) [(%V X ) - &oM] • (19) 

The renormalization factors Z m and Z are expected to be of order 0(1). The presence of 
the additive shift in the gaugino condensate bo(an) implies that the value of its jump at 
mug — is easier available than the value itself. 
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A first order phase transition should show up on small to moderately large lattices 
as metastability expressed by a two-peak structure in the distribution of some order 
parameter, in our case the value of the gaugino condensate. By tuning the bare parameters 
in the action, in our case the hopping parameter K for fixed gauge coupling /3, one can 
achieve that the two peaks are equal (in height or area). This is the definition of the phase 
transition point in finite volumes. By increasing the volume the tunneling between the 
two ground states becomes less and less probable and at some point practically impossible. 

In our simulations, besides the distribution of the gaugino condensate, we also studied 
other quantities as the string tension or the masses of the lightest bound states. The first 
results have been published recently fl3], together with a first hint for the existence 
of a phase transition from a simulation at (fl = 2.3, K = 0.195). In the present paper we 
keep the gauge coupling at (3 = 2.3 and exploit the region around K = 0.195. 

The Monte Carlo simulations are done by a two-step variant of the multi-bosonic 
algorithm ]T7| proposed in ||. We use polynomial approximations discussed in detail 
in [|TI| and correction procedures which are adapting some known methods from the 
literature fl9 , |20|1 to the present situation with Nf = | flavours. Our experience with 
this algorithm has been described already in previous publications [ |2"T| , |T^, |TB| and will 
be discussed in detail in a forthcoming paper p2|. 



Table 1: Parameters of the numerical simulations on 6 3 • 12 lattice at (3 = 2.3. 



K 


e 


A 


m 


n 2 


n 3 


n 4 


updates 


A 


Tplaq 


C (2A0) 


0.19 


0.0005 


3.6 


20 


112 


150 


400 


1487360 


0.888 


214(9) 


0.136(42) 


0.1925 


0.0001 


3.7 


22 


132 


180 


400 


3655680 


0.889 


220(7) 


0.220(36) 


0.195 


0.00001 


3.7 


24 


200 


300 


400 


460800 


0.892 


256(15) 


0.063(38) 


0.195* 


0.00003 


3.7 


22 


66 


102 


400 


1224000 


0.823 






0.196 


0.00001 


3.7 


24 


200 


300 


400 


952320 


0.889 


321(26) 


0.180(32) 


0.1975 


0.000001 


3.8 


30 


300 


400 


500 


506880 


0.926 


295(17) 


0.367(31) 


0.2 


0.000001 


3.9 


30 


300 


400 


500 


599040 


0.925 


317(16) 


0.424(26) 



The parameters of the numerical simulations on 6 3 • 12 lattice at j3 = 2.3 are sum- 
marized in table |l|. The run with an asterisk had periodic boundary conditions for the 
gaugino in the time direction, the rest antiperiodic. K is the hopping parameter and [e, A] 
is the interval of approximation for the first three polynomials of orders ^1,2,3, respectively. 
The fourth polynomial of order 714 is defined on [0, A]. In the eighth column the number of 
performed updating cycles is given. The ninth column contains the acceptance rate in the 
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noisy correction step A nc , the tenth column gives the exponential autocorrelation length 
for plaquettes r p i aq observed in the range of about 100 updating steps. The integrated au- 
tocorrelation is roughly a factor four higher, with large errors: for instance at K = 0.1925 
r piaq — 900 i 300. The last column contains the value of the autocorrelation function of 
the gaugino condensate at a distance 240, where the measurements were performed. 

The order parameter of the supersymmetry phase transition at zero gaugino mass is 
the value of the gaugino condensate 

P=fiE • (2°) 



X' 



The normalization is provided by the number of lattice points Q. We determined the 
value of p on a gauge configuration by stochastic estimators 

1 N v 

vEEfWVi (2i) 

ly V i=l xy 

on normalized Gaussian random vectors i] Xti . In practice N v = 25 works fine. Outside 
the phase transition region the observed distribution of p can be fitted well by a single 
Gaussian, but in the transition region a reasonably good fit can only be obtained with 
two Gaussians (see figure Jj). The fit parameters of the distributions (i = 1 or i = 1, 2) 

P- exp(-^#l (22) 



and the x 2 values per degrees of freedom are given in table 0. The normalization is such 
that P1+P2 = 1- Exact supersymmetry would imply that the widths of the two Gaussians 
are equal. This relation is broken by the lattice regularization and by the non-zero gaugino 
mass away from the phase transition point. In order to keep the number of fit parameters 
small we neglect this small symmetry breaking and the fits are done under the assumption 
or 1 = 02 = or. The statistical errors of the fit parameters are determined by jack-knifing 
64 statistically independent parallel runs. 

As figure [l] and table ^| show, in the region 0.195 < K < 0.1975 the distribution of the 
gaugino condensate can only be fitted well by two Gaussians. Comparing the two runs 
at K = 0.195 with antiperiodic, respectively, periodic boundary conditions in the time 
direction, one can see that the different boundary conditions do not have a sizeable effect 
on the distributions, as remarked before. For increasing K (decreasing bare gaugino mass) 
the weights shift from the Gaussian at larger p to the one with smaller p, as expected. 
The two Gaussians represent the contributions of the two phases on this lattice. The 
position of the phase transition on the 6 3 • 12 lattice is at K = 0.1955 ± 0.0005. The 
jump of the order parameter is Ap = p t — p 2 — 0.15. 
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Table 2: Fit parameters of the order parameter distributions corresponding to the 
runs in table The statistical errors in last digits are given in parentheses. 



K 


yi 


/ / 1 


1 In 




v 2 Id o f 

A /M.U.J. 


0.19 


1.0 


11.0023(26) 




0.0423(16) 


27.9/20 


0.1925 


1.0 


10.8807(30) 




0.0524(17) 


25.9/20 


0.195 


0.89(7) 


10.762(30) 


10.608(30) 


0.066(7) 


16.5/18 


0.195* 


0.83(6) 


10.78(3) 


10.60(3) 


0.055(7) 


16.3/18 


0.196 


0.35(7) 


10.722(11) 


10.588(11) 


0.073(3) 


5.7/18 


0.1975 


0.26(5) 


10.626(17) 


10.484(17) 


0.056(4) 


19.5/18 


0.2 


0.0 




10.3363(37) 


0.0562(18) 


21.4/20 



The two-phase structure can also be searched for in pure gauge field variables as the 
plaquette or longer Wilson loops. It turns out that the distributions of Wilson loops are 
rather insensitive. They can be well described by single Gaussians with almost constant 
variance in the whole range 0.19 < K < 0.2 (see, for instance, table ^j. This speaks 
against the appearance of a third chirally symmetric phase |23| , which has been suggested 
in \UA 



Table 3: Fit parameters of the plaquette distributions on 6 3 • 12 lattice at j3 = 2.3 
for different hopping parameters. The statistical errors in last digits are given in 
parentheses. 



K 


Pi 


Mi 


0~l 


X 2 /d.o.f. 


0.19 


0.974(25) 


0.63165(8) 


0.00425(13) 


0.89/47 


0.1925 


1.014(27) 


0.63511(8) 


0.00461(15) 


0.74/47 


0.195 


0.997(59) 


0.63811(19) 


0.00481(35) 


2.58/47 


0.196 


1.059(63) 


0.64182(22) 


0.00518(36) 


1.77/47 


0.1975 


0.987(54) 


0.64452(18) 


0.00444(30) 


2.26/47 


0.2 


1.018(44) 


0.64846(13) 


0.00424(22) 


2.00/47 
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4 Summary and discussion 



The observed dependence of the distribution of the gaugino condensate on the gaugino 
mass irtg near m- g = is consistent with a typical behaviour characteristic of a first order 
phase transition between two phases (see figure [I] and table 0). Our lattice volume (L 3 -T = 
6 3 • 12) is, however, still not very large in physical units, therefore the expected two-peak 
structure is not yet well developed. For instance, at K = 0.1925 we have LM®g ~ 3.6, 
with the smallest glueball mass M° [|13|, |2^|. In fact, a behaviour corresponding to a 



true first order phase transition can only be established in a detailed study of the volume 
dependence, which we postpone for future work. Therefore, the present observations are 
also consistent with a rapid cross-over at finite lattice spacings, approaching to a first 
order phase transition in the continuum limit (3 — > oo. On our 6 3 • 12 lattice for f3 = 2.3 
the phase transition (or cross-over) is at K = 0.1955 ± 0.0005. The jump of the gaugino 
condensate in lattice units is Ap ~ 0.15. 

A rather positive aspect of our Monte Carlo simulations is the ability of the two-step 
multi-bosonic algorithm || to cope with the difficult situation at small dynamical fermion 
mass in the environment of metastability of phases. 

In the numerical simulations we considered up to now only the unrenormalized gaugino 
mass and gaugino condensate. The transformation to the corresponding renormalized 
quantities defined in eqs. (|18| ) - (|I9|) will, however, not change the qualitative behaviour, 
because the multiplicative renormalization constants are expected to be of 0(1). One has 
to note that in the exploited range the bare gaugino masses m^g are small compared to the 
lightest bound state masses. With K = 0.1955 at K = 0.1925 we have m 0g /Mgg ~ 0.07. 
Similarly to QCD, it is expected that the mass gap in the spectrum is of the same order 
of magnitude as the scale parameter for the asymptotically free coupling A. As the 
preliminary results on the bound state masses show [|T3| , |22fl , at K = 0.1925 we already 
have an approximate degeneracy of the states which are expected to form the lowest chiral 
supermultiplet. 

Besides the volume dependence, another interesting question is the development of the 
phase transition signal towards the continuum limit at /3 — oo. In fact, the arguments in 
the introduction (at eqs. (0)-®) f° r the spontaneous chiral symmetry breaking Z 2 n c ~ ► Z 2 
refer to the continuum limit. The present numerical evidence shows that the discrete 
chiral symmetry breaking is manifested at non-zero lattice spacing in feasible numerical 
simulations and can be investigated by well established methods. 

Acknowledgements: It is a pleasure to thank Gemot Minister for helpful discus- 
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T3E-512 computer at HLRZ Julich. We thank HLRZ and the staff at ZAM for their kind 
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Figure 1: The probability distributions of the gaugino condensate for different hopping 
parameters at f3 = 2.3 on 6 3 • 12 lattice. The dased lines show the Gaussian components. 
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